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RESUMEN 

En este articulo se discuten los resultados de analisis inelasticos ante carga estatica monotona 
creciente que evaluan el comportamiento de edificios de 8 y 10 pisos estructurados con marcos 
ductiles momento resistentes. Los analisis no lineales se realizaron en OpenSees, donde la 
seccion transversal se modelo mediante fibras y elementos subdivididos en segmentos. Se 
consideraron siete patrones de carga lateral: tres relacionados con los modos de vibrar de cada 
edificio, dos de los modos acoplados, un patron rectangular y un patron de fuerzas inerciales 
que ocasiona la mayor deformacion en la azotea ante un acelerograma historico. 

PALABRAS CLAVE: marcos ductiles de acero, modo superior, curva de capacidad, distorsion 



INELASTIC RESPONSE OF STEEL BUILDINGS USING DIFFERENT 
LATERAL LOAD PATTERNS 



SUMMARY 

In this paper, results of pushover analyses are discussed in order to evaluate the behavior of 
buildings structured with moment-resisting ductile frames. Subject buildings are eight- and ten- 
stories in height. Nonlinear analyses were performed in OpenSees, where the cross-section was 
modeled with fibers and elements were sub-divided in segments. Seven lateral load patterns 
were considered: three related with the modes of vibration of each building, two related with the 
coupled modes, a rectangular pattern and a pattern associated with the inertial forces which 
develops the maximum roof deformation under an historical record. 
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1. INTRODUCCION 



Conforme a la tendencia reglamentaria internacional de diseno sismo resistente, las 
demandas se estiman considerando explfcitamente la respuesta inelastica de las estructuras; 
de manera que desde la etapa elastica se intenta predecir la respuesta plastica. Sin embargo, 
esto no es un proceso sencillo, puesto que esa prediccion conlleva simplificaciones y muchas 
variables intrinsecas dentro del proceso de analisis y diseno. 

La respuesta mas exacta y rigurosa del comportamiento inelastico que podria presentar 
una estructura se determina mediante analisis dinamicos no lineales de modelos complejos en 
programas de analisis especializados (Krawinkler 1998). Asf, la complejidad de la reproduccion 
de la respuesta dinamica de los elementos estructurales y el post-proceso dentro de los 
programas de analisis dinamicos no lineales, ocasionan que sea mucho mas comun que se 
prefiera realizar analisis no lineales ante carga estatica monotona creciente en la practica 
profesional y en la definicion de algunos lineamientos reglamentarios. Por ejemplo, la 
determinacion de los factores de reduccion para evaluar la respuesta de las edificaciones en el 
reglamento europeo EuroCode (Ap. 6.3.2 EC8-03, Tapia y Tena 2010). 

Conforme al FEMA 273 (NEHRP 1997), en el analisis no lineal ante carga estatica 
monotona creciente, un modelo que incorpora la respuesta inelastica del material se debe 
desplazar hasta un desplazamiento objetivo. Este modelo matematico del edificio se somete a 
un aumento monotono de fuerzas laterales o desplazamientos hasta que se supera el 
desplazamiento objetivo o se presenta el colapso de la estructura y pretende representar la 
demanda maxima de desplazamiento (y demandas) que podria experimentar el edificio durante 
un sismo real. Asf, siguiendo el analisis inelastico estatico se pretende que las fuerzas internas 
sean razonablemente aproximadas a las demandas que se presentarfan durante el sismo de 
diseno. Por esto, las demandas sfsmicas que se aplican al edificio en patrones laterales de 
carga deben estar relacionados con la distribucion de carga en la altura mas probable y mas 
crftica que se presentarfa durante un sismo intenso. 

Conforme al FEMA 273 (NEHRP 1997) para evaluar la capacidad inelastica de una 
estructura se deben considerar al menos dos patrones laterales de carga: 

a. Un patron uniforme basado en una proporcion de la masa de cada entrepiso 

b. Un patron modal. 

Este patron modal se define como un patron lateral que se determina mediante las 
cargas de entrepiso obtenidas por el metodo estatico (ecuacion 1) si la masa total participa en 
mas del 75 por ciento en el modo fundamental en la direccion en consideracion; este es el 
procedimiento mas comunmente utilizado. Sin embargo, el patron modal tambien puede ser 
definido como un patron lateral de carga proporcional a las fuerzas inerciales por entrepiso 
consistentes con la distribucion del cortante calculado por una combinacion modal usando: 



a. La respuesta del analisis del espectro del edificio incluyendo el suficiente numero de 
modos para considerar el 90 por ciento de la masa total. 

b. Un espectro apropiado del sitio. 



Fi =-W i h i 

1 Q 1 1 ZWihi 



( 1 ) 



En la ecuacion anterior, c es el coeficiente sfsmico de diseno, Q es el factor de 
comportamiento sfsmico, 1/17 es el peso de la Fesima masa y hi es la altura de la Fesima masa 
sobre el nivel de desplante (NTCS-04 2004). 
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En este articulo se evalua la respuesta inelastica ante diferentes patrones de cargas 
laterales de edificios de acero de mediana altura estructurados con marcos momento 
resistentes que fueron disenados para las mayores demandas sfsmicas que se propone en el 
Reglamento de Construcciones para el Distrito Federal Mexicano (RCDF-04) y que fueron 
proyectados siguiendo las constructivas mas empleadas en la Ciudad de Mexico. El proposito 
de la investigacion es evaluar los casos y las condiciones para las cuales los patrones de carga 
sugeridos por el FEMA 273 no permiten determinar la respuesta inelastica de una estructura. 



2. DESCRIPCION DE LOS MODELOS 

Los edificios estudiados, que tienen 8 y 10 pisos, estan estructurados con marcos 
momento resistentes de acero en una planta de 40 m por 35 m con cinco crujfas en cada 
direccion principal con anchos de 7.0 m y 8.0 m, respectivamente (figura 1). Se establecio que 
el edificio tendra un uso de oficinas para la determinacion de las cargas vivas y se supuso que 
estarfa ubicado en la subzona lllb conforme a la zonificacion sfsmica de la Ciudad de Mexico 
para asociarlo con las mayores demandas de aceleracion propuestas por el Reglamento de 
Construcciones para el Distrito Federal Mexicano (RCDF-04). 

Considerando las tendencias constructivas actuales, se supuso que las columnas 
estarfan construidas en secciones cajon con acero A572 Gr. 50 y que las trabes serfan con 
perfiles laminados IR con acero A572 Gr. 50; de manera, que el acero de ambos elementos 
tiene un esfuerzo de fluencia teorico de f y = 3,515 kg/cm 2 , un modulo de elasticidad igual a E= 
2,038,900 kg/cm 2 . El sistema de piso supuesto es lamina acanalada tipo Losacero; sin 
embargo, su contribucion como diafragma se desprecio en este estudio. 

Las cargas muertas y cargas vivas se determinaron mediante las Normas Tecnicas 
Complementarias sobre Criterios y Acciones (NTCCr-04) del Reglamento de Construcciones 
para el Distrito Federal Mexicano (tabla 1); mientras que la determinacion de la seccion 
transversal de los elementos estructurales se realizo mediante un proceso iterativo de diseno en 
el programa ETABS (CSI 2010). Es posible consultar mayor informacion sobre las cargas 
actuantes consideradas y los pormenores del proceso de diseno por capacidad en Gutierrez 
(2013). 




Figura 1. Configuracion geometrica de los edificios. 
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Tabla 1. Cargas distribuidas consideradas en este estudio 



Solicitacion 


Entrepiso 

(MPa) 


Azotea 

(MPa) 


Carga muerta 


6.57 


5.69 


Carga viva 


2.45 


9.81 


Carqa viva accidental 


1.77 


6.87 



El espectro de diseno se determino considerando el criterio del Apendice A de las 
Normas Tecnicas Complementarias para el Diseno por Sismo (NTCS-04) del Reglamento de 
Construcciones para el Distrito Federal Mexicano (RCDF-04), donde se supuso un periodo 
dominante del suelo donde se localizarfa la estructura es igual a T s = 2.0 segundos en 
congruencia a las condiciones de la zona sismica lllb (zona del ex-Lago con arcillas altamente 
compresibles). El diseno de los edificios se realizo considerando la maxima ductilidad permitida 
en el RCDF-04 para este sistema estructural mediante la aplicacion de un factor de 
comportamiento sismico igual a Q= 4. El factor de reduccion por sobrerresistencia se determino 
conforme al criterio propuesto en el Apendice A (ecuacion 2, figura 2a). Es posible consultar una 
discusion mas amplia en Tapia y Tena (201 0). 



R = 



10 






SiT <T a 
SiT >T a 



( 2 ) 



En la ecuacion anterior, T es el periodo de la estructura y T a es el periodo caracterfstico 
del espectro de diseno que en el RCDF depende del periodo dominante del suelo T s (NTCS-04). 
Asi, en la figura 2b se muestra el espectro de diseno elastico y el espectro inelastico 
considerado en esta investigacion. 





Periodo (seg) 



a) Factor de reduccion por sobrerresistencia b) Espectros de diseno (zona lllb) 

Figura 2. Consideraciones para el diseno sismico segun NTCS-04 



Siguiendo las tendencias de la practica en Mexico, en el estudio se realizo un cambio 
de seccion transversal entre las columnas perimetrales del edificio y las columnas de los 
marcos interiores. En la tabla 2 se muestran las secciones estructurales finales consideradas en 
este estudio despues del proceso iterativo de diseno. 

Ademas, en la figura 3 se muestra las relaciones de esfuerzo por entrepiso (demanda 
entre resistencia) de las trabes y las columnas (ecuacion 3) para hacer notar que en el diseno 
se observo que las trabes constituyeran la primera Ifnea de defensa en la respuesta inelastica y 
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que las columnas se disenaron para que tuvieran poco o nulo dano (mecanismo de colapso viga 
debil - columna fuerte). 



Tabla 2. Secciones transversales finales consideradas en este estudio 



Elemento 

estructural 


Modelo de 8 pisos 


Modelo de 10 pisos 


Entrepiso 


Seccion 


Entrepiso 


Seccion 


Trabes 


N1 a N4 
N5 a N7 
N8 


IR 27"x140 kg/m 
IR 24"x140.1 kg/m 
IR 24"x82 kg/m 


N1 a N5 
N6 a N9 
N10 


IR 24"x125.1 kg/m 
IR 24"x140.3 kg/m 
IR 1 8"x52.2 kg/m 


Columnas de 
marcos interiores 


N1 a N3 
N4 a N6 
N7 a N8 


OR 60x60 cm; t= 3.18 cm 
OR 50x50 cm; t= 2.85 cm 
OR 50x50 cm; t= 1 .91 cm 


N1 a N3 
N4 a N7 
N8 a N10 


OR 60x60 cm; t=3.49 cm 
OR 60x60 cm; t= 2.85 cm 
OR 50x50 cm; t= 2.54 cm 


Columnas de 

marcos 

perimetrales 


N1 a N3 
N4 a N6 
N7 a N8 


OR 60x60 cm; t= 2.85 cm 
OR 50x50 cm; t= 2.22 cm 
OR 50x50 cm; t= 1 .91 cm 


N1 a N3 
N4 a N7 
N8 a N10 


OR 60x60 cm; t=3.18 cm 
OR 60x60 cm; t= 2.54 cm 
OR 50x50 cm; t= 1 .91 cm 



e = 



p u _|_ M ux _|_ M uy < ^ Q 
Pr M Rx M Ry 



( 3 ) 



En la ecuacion 3, que es aplicable para miembros en flexo-compresion (IMCA 2014), P u 
es la carga axial actuante; M ux y M uy son los momentos actuantes en ambas direcciones; Pr es 
la carga axial resistente y Mr x y MRy son los momentos resistentes en las direcciones 
principales. La ecuacion debe cumplir que e < 1.0 para asegurar que la seccion transversal 
propuesta satisface las solicitaciones por flexo-compresion. 



13 8 
> 

Z 7 

6 

5 



-OO 003K»— <X»- 



A AflMPOOO 



cesxo 



A-i 

-A- 



00®0 <3XCOO«> 



»o oo 



<x><»o cco 



4 -- A 4YS*W£>0 



3 -- A AA3CS>ASI* 
2 A - <X> 



I 

_±_ 

I 



<X> CSXHSOO 



OO <2>OCCSX> 



L AK^I — CS>- 



4 t OO <KXO«X> 

I I 



_i i i i_ 



0,00 
O Vigas 
A Columnas 



0,25 



0,50 



0,75 



1,00 



a) Modelo de 8 pisos 
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b) Modelo de 10 pisos 



Figura 3. Relaciones de esfuerzo (demanda contra resistencia) de los modelos 
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3. MODELADO DE ANALISIS INELASTICO 

Se realizaron analisis inelasticos ante carga estatica monotona creciente a los modelos 
en el programa de analisis OpenSees, Open System for Earthquake Engineering Simulation 
(Mazzoni et at 2006). Este programa es un software que usa el metodo de elemento finito para 
desarrollar aplicaciones que simulan el comportamiento estructural y geotecnico de sistemas 
sujetos a excitaciones sfsmicas, desarrollado por el Centro de Investigacion de Ingenierfa 
Sfsmica del Pacifico (Pacific Earthquake Engineering Research Center , PEER). 

En el modelo de OpenSees, se uso el material Stee/02 (el modelo de Giuffre- 
Menegotto-Pinto) con endurecimiento por deformacion isotropico en los perfiles de acero. 
Ademas, debido a la respuesta inelastica del material se propuso una solucion iterativa y que se 
resuelve con el algoritmo de Newton a partir de una librerfa del programa. Las ecuaciones se 
formaron mediante un esquema general, que se numero usando el metodo Cuthill-Mckee, que 
optimiza la numeracion de los nodos para reducir el ancho de banda de la matriz. Este 
procedimiento despliega un mensaje cuando la estructura se desconecta (colapsa). 

Con el proposito de considerar los efectos de segundo orden, los elementos trabes y 
columnas se subdividieron en ocho sub segmentos como se muestra en la figura 4a. Ademas, 
las columnas se modelaron considerando un desplome originado por causa de su fabricacion 
y/o montaje, mediante una deformacion al centro de su longitud; como se muestra en la figura 
4b. La deformacion relacionada con la imperfeccion se considero como la maxima magnitud que 
el Instituto Mexicano de la Construccion en Acero (IMCA 2014) sugiere como permisible 
(ecuacion 4). 



Desplome m £ X — — — 0.001L 



(4) 





• i 

<-► 

0.001L 



a) Discretizacion de los elementos en segmentos 



b) Imperfeccion en 
columnas 



Figura 4 Detalles del modelo inelastico en OpenSees 

En la figura 5a se muestra la configuracion tipica del modelo estudiado que tiene el 
proposito de estimar la rigidez lateral de los edificios. Los marcos del edificio se unieron 
mediante elementos de alta rigidez que se modelaron elasticos y cuya unica funcion es 
transmitir la carga lateral. 
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Columnas 



Trabes 

b) Seccion transversal 



tttt 



t-V 



ttti 



Figura 5. Consideraciones de los modelos inelasticos 



Con el proposito de considerar la distribucion de la plasticidad en la seccion transversal 
de los elementos (figura 5b), las secciones transversales en el modelo se propusieron con fibras 
de longitud finita con ocho puntos de integracion por elemento para modelar la variacion de la 
deformacion a lo largo de cada miembro. La seccion transversal de las vigas de seccion I, se 
dividio en tres parches, dos para los patines y una para el alma, los cuales se discretizaron con 
fibras con forma cuadrilatera. Mientras que las columnas de seccion cajon se discretizaron 
mediante cuatro fibras rectangulares para formar las paredes laterales de la seccion cajon con 
seis subdivisiones horizontales y cinco verticales y cuatro secciones mas para las esquinas de 
la seccion cajon. A fin de evitar inestabilidad por torsion, las fibras se unieron con el comando 
section Aggregator que aporta la rigidez por torsion a las secciones subdivididas en fibras. 



4. PROTOCOLOS DE CARGA 

En los analisis inelasticos ante carga estatica monotona creciente, se determinaron siete 
patrones de carga lateral como se describe en los siguientes apartados: un patron lateral 
rectangular, tres relacionados con los modos de vibrar de cada estructura, dos patrones que 
estan asociados al acoplamiento de los tres primeros modos y un patron relacionado con las 
fuerzas inerciales que provocan la mayor deformacion en la azotea en un analisis dinamico ante 
el acelerograma historico registrado en durante el sismo de septiembre de 1985 en la Ciudad de 
Mexico. 



4.1 Patron uniforme 

Con base en las recomendaciones del FEMA 273 (NEHRP 1997) se propuso un patron 
lateral de cargas constantes en la altura; es decir, un patron rectangular uniforme. La magnitud 
de las cargas se normalizaron con la fuerza lateral en la azotea de los edificios obtenidas del 
metodo estatico de las NTCS-04 (ecuacion 1). 
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4.2 Patrones modales 

Se determinaron patrones de carga relacionados con los tres primeros modos de vibrar 
siguiendo el procedimiento que se ejemplifica a continuacion para el modelo de 10 niveles. En 
los analisis, la matriz de rigidez lateral de los edificios se obtuvieron utilizando las formulas 
Wilbur. Es posible encontrar mas informacion sobre los criterios para determinar la rigidez 
lateral, resistencia y balance de cargas en estructuras de acero en Gutierrez (2013) y Tapia et 
at. (2013). 

La rigidez y el desplazamiento del primer entrepiso de una estructura cuando se supone 
que las columnas empotradas en la cimentacion esta definida por la ecuacion 5; mientras que la 
rigidez y el desplazamiento del segundo piso considerando las mismas condiciones de apoyo 
estan definidos por la ecuacion 6. 



Ri 



48 E 
(Olftl) 



; d 1 



( i N (h 1 + h 2 ) 

VzkJ (EK tl #§a) 



( 5 ) 



| = 48J? . n _ ( 4fa 2 \ , (fti +h 2 ) Ch 2 +h 3 ) 

2 [p 2 h 2 y U2 kxkcJ + (j iKtl+ U££iy Y.K t2 



( 6 ) 



En las ecuaciones, E es el modulo de elasticidad, Rn es la rigidez del entrepiso n, Ktn es 
la rigidez de las vigas del nivel sobre el entrepiso n\ K cn es la rigidez de las columnas del nivel 
sobre el entrepiso n y h n es la altura del entrepiso n. En las ecuaciones, la rigidez de la trabe 
(K tn = It/U) y la rigidez de las columnas (K cn = lc/h c ) se calculan con el momento de inercia de las 
trabes It y el momento de inercia de las columnas l c , en funcion de la longitud neta de las trabe 
L t y la altura neta de las columna h c . Es posible encontrar mas informacion sobre estas 
ecuaciones en Tapia et ai (2013). 

Las condiciones de frontera del apoyo de las columnas del edificio no tiene influencia 
sobre los entrepisos superiores, por lo que la rigidez y el desplazamiento en los siguientes 
entrepisos pueden ser determinados aplicando la ecuacion 7. En ella, m, n y o son subindices 
que identifican los tres pisos consecutivos de abajo hacia arriba. 



48£~ _ jj f 4h n \ , th-m+hn) , n + h 0 ) 

(Dnh-n) n ^ L l<ciJ (L Ktm) L Ktn 



Mediante este procedimiento y considerando las secciones transversales de la tabla 2, 
fue posible determinar las magnitudes de las rigideces de los entrepisos. En la tabla 3 se 
muestran las magnitudes obtenidas para el modelo de diez pisos, mientras que la matriz de 
rigidez, la matriz de masas [M] y la matriz de rigidez lateral [K]i_at se muestran en las ecuaciones 
8, 9 y 10, respectivamente para ejemplificar el procedimiento. 

Los modos de vibrar de las estructuras estudiadas se obtuvieron mediante la resolucion 
de un problema de valores unicos. La ecuacion de movimiento para un sistema de n grados de 
libertad esta definida por la ecuacion 1 1 cuando el sistema no es amortiguado. 
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Tabla 3. Rigideces laterales del modelo de 1 0 niveles 



Entrepiso 


Rigidez K n (kg cm) 


10 


kio 


3.07 (10) 5 


9 


kg 


4.81 (10) 5 


8 


k8 


4.81 (10) 5 


7 


k 7 


5.21 (1 0) 5 


6 


k6 


5.52 (10) 5 


5 


ks 


5.87 (10) 5 


4 


k 4 


5.87 (10) 5 


3 


k3 


5.96 (1 0) 5 


2 


k 2 


7.24 (10) 5 


1 


ki 


1.27 (10) 6 



\( k i 



' Lat 



~ k 2 { k 2 + k 3) ~ k 3 

~ k 3 ( k 3 + k 4 ) 



~ k 4 

0 

0 

0 

0 

0 

0 



{ k 4 + k 5 



{ k 5 



[ k 6 + k 7 ) 



[ky + kg') 



- k & 

0 

0 



0 

0 

0 

0 

0 

0 

—ko 



( k 8 + k 9 ) 



~k n 



0 

0 

0 

0 

0 

0 

0 

—kn 



(kg + k 



10) 



*10 



0 

0 

0 

0 

0 

0 

0 

0 

~ k 10 

k 10 



( 8 ) 



1113.15 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1100.13 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1100.13 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1086.21 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1086.21 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1081.73 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1081.73 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1062.05 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


1062.05 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


908.31 



1.99E+06 


-7.24E+05 


0 


0 


0 


0 


0 


0 


0 


0 


-7.24E+05 


1.32E+06 


-5.96E+05 


0 


0 


0 


0 


0 


0 


0 


0 


-5.96E+05 


1.18E+06 


-5.87E+05 


0 


0 


0 


0 


0 


0 


0 


0 


-5.87E+05 


1.17E+06 


-5.87E+05 


0 


0 


0 


0 


0 


0 


0 


0 


-5.87E+05 


1.14E+06 


-5.52E+05 


0 


0 


0 


0 


0 


0 


0 


0 


-5.52E+05 


1.07E+06 


-5.21E+05 


0 


0 


0 


0 


0 


0 


0 


0 


-5.21E+05 


1.00E+06 


-4.81E+05 


0 


0 


0 


0 


0 


0 


0 


0 


-4.81E+05 


9.61E+05 


-4.81E+05 


0 


0 


0 


0 


0 


0 


0 


0 


-4.81E+05 


7.88E+05 


-3.07E+05 


0 


0 


0 


0 


0 


0 


0 


0 


-3.07E+05 


3.07E+05 



[M]{x}+[K]{x} = { 0} 



( 11 ) 
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Donde, [M]es la matriz de masas, [K]i_at es la matriz de rigidez lateral, {x} es el vector de 
aceleraciones y {x} es el vector de desplazamientos. Ahora, con el proposito de obtener los 
modos de vibracion, es necesario resolver un sistema de ecuaciones diferenciales lineales 
ordinarias, por lo que su solucion se obtiene mediante una funcion exponencial y utilizando la 
transformada de Laplace. Asi, el vector de desplazamientos queda definido por la ecuacion 1 2 y 
es posible demostrar que el vector de aceleraciones esta definido por la ecuacion 13; donde {A} 
es un vector de amplitudes. 



*2 




*1 

*2 




( 12 ) 



(13) 



Asi, sustituyendo los vectores anteriores, considerando el cambio de variable {-X= s 2 ) y 
despreciando la solucion trivial ({A}={0}), es posible demostrar que la solucion para un sistema 
de varios grados de libertad esta definida como se indica en la ecuacion 14. 



\\K uu -m\ = 0 ( 14 ) 

En esta investigacion se utilizo el comando “genvecs (K, M)” del programa Mathcad 14 
(PTC 2007) que resuelve los calculos matematicos para determinar los eigen valores 
relacionados con los modos de vibracion. La configuracion de los primeros tres modos de vibrar 
y la magnitud de los periodos se ilustran en la figura 6. 
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a) Modelo de ocho niveles b) Modelo de diez niveles 

Figura 6. Primeros modos de vibrar de los modelos estudiados 
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4.3 Patrones de analisis dinamicos 

En este estudio se considero un patron de carga que estuviera relacionado con las 
fuerzas inerciales por entrepiso generadas por un sismo historico; para lo cual se realizaron 
analisis dinamicos no lineales a los modelos de ambos edificios bajo la excitacion provocada 
por el acelerograma registrado en la estacion de la Secretaria de Comunicaciones y 
Transposes en la Ciudad de Mexico durante el sismo de septiembre de 1985 en la direccion 
este - oeste (SCT). 

En la figura 7a se muestra el espectro de pseudo-aceleracion del registro seleccionado 
que fue post-procesado en el programa NonLin (Charney 2011); mientras que en la figura 7b se 
muestra la historia del desplazamiento en la azotea cuando el modelo de 10 pisos es sometido 
al acelerograma seleccionado. En los analisis, la maxima deformacion en la azotea ocurrio a los 
56.46 seg y 57.88 seg en los modelos de ocho y diez niveles, respectivamente. Entonces, se 
determinaron los cortantes resistentes por entrepiso en las columnas asociados con ese paso 
para definir los patrones de carga lateral que se discuten mas adelante. 
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Figura 7. Determinacion del patron de fuerzas inerciales del modelo de 1 0 pisos 



Asf, en los analisis inelasticos paso a paso de ambos modelos se determinaron las 
fuerzas inerciales (cortante lateral por entrepiso Vi) relacionado con esa demanda maxima. 
Finalmente, el patron de carga lateral se obtuvo normalizando el cortante lateral en la azotea 
producido por el registro historico con la fuerza lateral obtenida del metodo estatico. 



4.4 Patrones de los modos acoplados 

Una vez que se determinaron los modos de vibrar, se calcularon patrones acoplados 
mediante una combinacion de los modos usando el metodo de la rafz cuadrada de la suma de 
los cuadrados (metodo SRSS por sus siglas en ingles) siguiendo la ecuacion 15. 

R = ( 15 ) 

En la ecuacion n es el numero de modos considerados, / es el numero de modo y Ft la 
respuesta combinada. Siguiendo este procedimiento, se obtuvieron dos patrones: un patron 
considerando solo los modos 1 y 2 que se titulo M1+M2 y otro patron que se calculo teniendo en 
cuenta los modos 1 , 2 y 3 que se titulo M1+M2+M3. 
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En el procedimiento descrito se combinan cargas laterales aleatorias sin incluir 
explicitamente la contribution de la masa modal efectiva de cada modo. Es posible encontrar 
formulas simplificadas para incluir esos efectos en Christopoulos y Filiatrault (2006). 



5. CORTANTE DE ANALISIS DINAMICO 

Los patrones de carga considerados en este estudio se normalizaron a partir de la 
magnitud del patron de cargas laterales del metodo estatico convencional de las NTCS-04 en el 
nivel de la azotea (ecuacion 1). En la figura 8 se muestra la configuration de los patrones de 
carga normalizados con exception del modelo rectangular que es constante en toda la altura. 
Estos patrones de carga lateral concuerdan con los resultados de otras investigaciones como 
Chopra y Goel (2002) y Chopra et at. (2004). 







Patron modo 2 



Patron modo 3 Patron SCT 
a) Modelo de ocho pisos 



Patron M 1 +M 2 






M1+M2+M3 



0,0 0,5 1,0 




Patron modo 2 



Patron M 1 +M 2 



Patron modo 3 Patron SCT 
b) Modelo de diez pisos 

Figura 8. Patrones de carga normalizados considerados en esta investigation 



M1+M2+M3 



Particularmente, el patron de distribucion SCT con forma triangular parecerfa tener 
relation con la distribucion del metodo estatico, pero no es asf; se trata de una distribucion 
relacionada con una demanda dada. 
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6. ANALISIS INELASTICOS 

En la figura 9 se muestra la configuracion deformada esquematica del modelo de ocho 
pisos en los analisis inelasticos estaticos ante carga monotona creciente con el patron de 
cargas relacionado con el modo fundamental para hacer notar que la configuracion deformada 
de los elementos y rotaciones en los nodos de conexion tiene congruencia con la respuesta 
esperada de un marco momento - resistente. 




Con el proposito de estudiar la evolucion del dano de los patrones de carga 
relacionados con la configuracion modal, en la figura 10 se muestra el comportamiento de las 
distorsiones de entrepiso en el modelo de 10 niveles en el paso donde la distorsion global en el 
analisis ante las cargas laterales asociadas al modo principal (modo 1) tiene una magnitud igual 
a 0.01 H, 0.02Hy 0.03H. En otras palabras, en la figura se muestra la configuracion deformada 
relacionada con la magnitud del patron de carga (el paso) que fue necesario aplicar para que en 
el analisis ante el patron de cargas del modo fundamental se presentara esa magnitud de 
deformacion en la azotea. 




Modo 1 Distorsion (%) 

• Modo 2 
Modo 3 

c) Distorsion global 0.01H 




Modo 1 Distorsion (%) 

• Modo 2 
Modo 3 

b) Distorsion global 0.02H 




Modo 1 Distorsion (%) 

• Modo 2 
• • • • • Modo 3 

c) Distorsion global 0.03H 



Figura 10. Evolucion de las deformaciones en el modelo de 10 pisos 



Asf por ejemplo, la distorsion del modo 1 cuando la distorsion global es igual a 0.01H 
ocurre cuando el patron de cargas es afectado por un factor igual a 0.46; de manera que en la 
comparacion se uso la deformacion que se presentaba en los analisis ante los patrones de 
carga relacionados con el modo 2 y modo 3 cuando el factor incremental era proximo a 0.46. 
Siguiendo este procedimiento en cada caso, como se ilustra en la tabla 4, se obtuvieron las 
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distorsiones de entrepiso que se muestran en las figura 9. En estos calculos, la distorsion global 
se determino como la deformacion en la azotea dividida entre la altura total del modelo H. 

Se nota el patron de carga relacionado con el primer modo lidera las deformaciones 
desde el inicio del analisis (figura 9). Esto implica que el patron de cargas del metodo estatico 
convencional de las NTCS-04, esta relacionado igualmente con las mayores demandas de 
deformacion; ya que ese criterio (ecuacion 1) se obtiene considerando una distribucion 
triangular inversa que corresponde al modo 1. Esta observacion concuerda con otros resultados 
(FEMA 440 2005) de analisis de varios modelos con distintos patrones de cargas y analisis 
modales. 

Segun el Apendice A de las NTCS-04, las distorsion permisible para la revision del 
estado limite de seguridad contra colapso es igual a 0.03 h, donde h es la altura del entrepiso. 
Con base en estos resultados, los modelos de este estudio son capaces de desarrollar una 
magnitud superior a la permisible sin colapsar. Esta observacion coincide con otros estudios 
como Tapia y Tena (2013), donde modelos de edificios racionalmente disenados siguiendo 
estrictamente la metodologfa de diseno de las NTCS son capaces de desarrollar deformaciones 
(y ductilidades) superiores a las establecidas en los reglamentos especializados vigentes. 



Tabla 4. Comparacion de la deformacion de entrepiso en distintos pasos (modelo de 1 0 pisos) 



Patron de carga 


Modo 1 


Modo 2 


Modo 3 


Deformacion 


0.01H 


0.02H 


0.03H 


0.01H 


0.02H 


0.03H 


0.01H 


0.02H 


0.03H 


Factor de 
incremento 


0.46 


0.91 


1.22 


0.46 


0.92 


1.23 


0.45 


0.91 


1.21 


10 


36.12 


72.12 


108.12 


26.12 


52.12 


70.12 


24.12 


48.12 


64.12 


9 


31.41 


62.69 


93.86 


22.48 


44.88 


60.36 


20.81 


41.53 


55.33 


8 


26.83 


53.51 


79.99 


19.08 


38.09 


51.21 


17.84 


35.62 


47.45 


7 


22.30 


44.45 


66.36 


15.94 


31.83 


42.79 


14.98 


29.93 


39.88 


6 


17.89 


35.63 


53.17 


13.03 


26.03 


35.00 


12.13 


24.23 


32.29 


5 


13.63 


27.12 


40.52 


10.22 


20.42 


27.48 


9.34 


18.66 


24.88 


4 


9.66 


19.21 


28.79 


7.49 


14.97 


20.19 


6.78 


13.56 


18.10 


3 


6.13 


12.21 


18.45 


4.92 


9.84 


13.32 


4.48 


8.96 


11.99 


2 


3.23 


6.45 


9.90 


2.67 


5.35 


7.30 


2.46 


4.93 


6.63 


1 


1.11 


2.25 


3.55 


0.95 


1.91 


2.67 


0.89 


1.79 


2.45 



En la figura 1 1 se muestran las curvas de capacidad obtenidas en los analisis 
inelasticos ante carga estatica monotona creciente mediante los siete patrones de carga lateral 
considerados en esta investigacion. Las curvas muestran el comportamiento del cortante basal 
calculado como la demanda sobre las columnas contra la distorsion global calculado con la 
deformacion en la azotea dividida entre la altura total del modelo. Las graficas muestran una 
disminucion muy gradual de la rigidez lateral, que es un comportamiento tfpico de modelos que 
consideran una plastificacion detallada mediante fibras y elementos subdivididos en varios 
segmentos (Tapia y Tena 2013). 

El patron de cargas asociado al modo fundamental rige el diseno (figura 11), debido a 
que se requiere una menor magnitud de la carga lateral para lograr que la estructura tenga un 
mayor desplazamiento. En el modelo de 8 pisos, este patron (modo 1) estima adecuadamente 
las demandas maximas que se podrfan presentar en un sismo, representadas por la curva del 
patron de las fuerzas inerciales (patron SCT). 
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Sin embargo, ante el aumento de la altura del edificio (modelo de 10 pisos), las 
demandas relacionadas con el patron del modo 1 sobreestiman las demandas maximas que 
reporta la demanda del sismo historico. La curva de capacidad del modo acoplado M 1 +M 2 
tiende a los resultados del patron del modo 1 y denota, igualmente, una dependencia con la 
altura de los edificios. 




M °d° 1 Distorsion global (%) 

— — Modo 2 
• • • • Modo 3 

‘Rectangular 

M1+M2 

M1+M2+M3 




Modo 1 Distorsion (%) 

Modo 2 

• • • • • Modo 3 

‘Rectangular 

M1+M2 

M1+M2+M3 



a) Modelo de 8 niveles b) Modelo de 10 niveles 

Figura 11. Curva de capacidad global de los modelos estudiados 



En general, el comportamiento ante los patrones del modo 2, modo 3, el modo acoplado 
M1+M2+M3 y el patron rectangular son muy similares y sus demandas son poco significativos en 
los edificios estudiados (figura 11). Esto significa que en los casos estudiados, los modos 
superiores no tienen una participation importante en la respuesta inelastica. 



7. CONCLUSIONES 

Esta investigacion estudia el comportamiento inelastico de dos edificios de ocho y diez 
niveles con uso de oficinas estructurados con marcos momento resistentes de acero, que se 
supusieron localizados en la zona lllb de la subzonificacion sfsmica del Distrito Federal 
Mexicano (Zona del ex-Lago) para asociarse a las mayores demandas de aceleracion 
propuestas en las Normas Tecnica Complementarias para el Diseno por Sismo (NTCS-04) del 
Reglamento de Construcciones para el Distrito Federal (RCDF-04) y se disenaron para la 
maxima condition de ductilidad permitida para este sistema estructural con un factor de 
comportamiento sfsmico Q= 4.0. 

Esta investigacion pretende evaluar si los patrones de carga lateral que propone el 
FEMA 273 y que son ampliamente aceptados en la ingenieria son una representation 
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razonablemente aproximada de las demandas que se presentarian durante el sismo de diseno 
considerando siete patrones laterales: 

a. Tres patrones modales asociados a la configuracion del modo fundamental (un patron 
uniforme basado en una proporcion de la masa de cada entrepiso), del modo 2 y del 
modo 3. 

b. Un patron rectangular. 

c. Un patron lateral de fuerzas inerciales de entrepiso (cortante lateral) relacionados con la 
maxima deformacion en la azotea de los modelos ante el acelerograma registrado en 
SCT en el sismo de septiembre de 1985. 

d. Dos patrones relacionados con el acoplamiento de los modos: uno de ellos 
considerando el primer y segundo modos y otro considerando los tres primeros modos. 

Los analisis no lineales se realizaron en OpenSees mediante un modelo refinado que 
describe la plastificacion mediante una discretizacion de las secciones transversales con fibras 
y una subdivision de los elementos en segmentos para incluir detalladamente los efectos de 
segundo orden. Ademas, las columnas se modelaron con una deformacion inicial U1000 al 
centro del claro que se establecio con la maxima tolerancia por errores de montaje que 
establecen los reglamentos vigentes especializados en estructuras de acero. 

De acuerdo con los resultados de la investigacion, a partir de la evolucion de las 
distorsiones y de curvas de capacidad se pueden hacer las siguientes observaciones: 

• La evolucion de las distorsiones de los entrepisos ante los patrones de carga 
relacionados con la configuracion deformada de los modos 1, 2 y 3 demuestran que, 
desde el inicio de los analisis, el modo fundamental representa la mayor demanda para 
una misma magnitud del factor incremental del patron de carga. 

• Los resultados de los analisis no lineales ante los patrones laterales considerados estan 
relacionados con deformaciones grandes sin que se presente el colapso. Esto 
demuestra que las estructuras de acero son flexibles y, en todos los casos, superan la 
distorsion permisible para la revision del estado limite de colapso propuesta en las 
NTCS-04 para marcos ductiles de acero (3%). 

• Con base en las curvas de capacidad, el patron de cargas asociado al modo 
fundamental rige el diseno, debido a que se requiere una menor magnitud de la carga 
lateral para lograr que la estructura tenga un mayor desplazamiento. 

• En el modelo de 8 pisos, el patron del modo fundamental es una buena representacion 

de las demandas que generarfa un sismo, representado por el patron de fuerzas 

inerciales de un registro historico (SCT). Sin embargo, ante el aumento de la altura en 
el modelo de 10 niveles, el patron del modo fundamental sobreestima las demandas. 

• Las demandas que producen los patrones del modo 2, modo 3, el modo acoplado 

M1+M2+M3 y el patron rectangular son similares y requieren una magnitud mayor para 

danar las estructuras estudiadas, por lo que el patron del modo acoplado M1+M2 es mas 
crftico. Esto significa que los modos superiores tienen poca influencia en la respuesta 
inelastica de los edificios estudiados. 

Los autores hacen notar que aun es necesario ampliar el estudio para incluir un mayor numero 
de modelos antes de generalizar los resultados obtenidos. Pese a que en este estudio se 
considero la tendencia de la practica y a que los modelos no son perfectamente regulares en 
planta y se incluyeron deformaciones iniciales, es posible que otros modelos esten relacionados 
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con resultados significativamente diferentes. Los resultados dependen altamente de la forma de 
la edificacion, sistema estructural, consideraciones del modelado, conexiones, distribucion de 
masas y sistemas resistentes, entre otros. Se sugiere realizar estudios adicionales para 
establecer los casos, Ifmites y condiciones que correspondan a los resultados de este artfculo. 
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